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C/5 

We discuss scaling laws of fusion yields generated by laser-plasma interactions. The 
yields are found to scale as a function of the laser power. The origin of the scaling 
law in the laser driven fusion yield is derived in terms of hydrodynamic scaling. We 
point out that the scaling properties can be attributed to the laser power dependence of 
O three terms: the reaction rate, the density of the plasma and the projected range of the 

* ^ plasma particle in the target medium. The resulting scaling relations have a predictive 

power that enables estimating the fusion yield for a nuclear reaction which has not been 
investigated by means of the laser accelerated ion beams. 

_ 1. Introduction 

cn 

>>■ Ultraintense laser pulses can generate high energy protons when irradiating a thin 

foil target 1234^ r^^ie proton beams generated by the laser irradiations have var- 
ious potential applications - in science, engineering and medical imaging m The 
10 laser pulse parameters, i.e., pulse energy, peak intensity and pulse duration, are, 

^\ however, currently not optimized for practical applications. The determination of 

scaling laws in both the characteristics of the accelerated ion beams itself and in the 
' I resulting fusion yield is crucial for such an optimization. In this respect, we restrict 

our study on the scaling laws of fusion yields derived empirically from the currently 
available experimental data. We do not consider the acceleration mechanism of re- 
diation pressure acceleration 1^ that could be important to obtain higher energy ion 
^ beams for the oncological proton(or hadron)-therapy f^, because such a mechanism 

is not feasible at current technology. Clarifying the scaling laws of the fusion yield 
promises to extend the result even to estimate fusion yield for not well-investigated 
nuclear reactions by means of the ultrahigh-intensity laser-matter interactions, e. 
g., aneutronic fusion reactions IMHI and deuteron-induced reactions 13I1IMIH Es- 
pecially deuteron-induced reactions are expected to induce enough yields for the 
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radioactive isotope production for PET diagnostics using laser-induced ions. For 
the purpose of the radioactive isotope production for PET diagnostics, monochro- 
matic ion beam with an energy as high as 200 MeV for oncological hadron-therapy 
is not necessary. 

On this issue, scaling laws of the maximum energy and conversion efficiency of 
laser-accelerated protons from thin foil targets have been studied | i5 | i6 | i7 [ -yyi^en 
a laser of intensity higher than 10^^ W/cm^ is irradiated on an aluminum foil 
target, a simple power scaling of the form Epmax = a x l'' with b = 0.5 ± 0.1 is 
found null, where Epmax and / denote the maximum proton energy and the laser 
pulse intensity, respectively. More in general the maximum proton energy is derived 
also as a function of the laser power ^^3. The dependence of the conversion efficiency 
on the laser energy is found to be linear, as it is shown in Fig. 2 in 

The scaling feature of the accelerated protons is obviously reflected in the yield 
of the laser-induced nuclear reactions. The neutron yield in deuterium cluster ex- 
plosions is found to scale as a function of the laser pulse energy (E^) IHIIll and it 
follows a quadratic power law dependence. The other examples are the ^^Zn yield 
and the ^^C yield^^Slfor the reactions ^^Cu(p, n)^^Zn and ^^Ti(ji, n)^^C, respectively. 
Both reactions are obtained with proton beams generated at the VULCAN laser of 
the Rutherford Appleton Laboratory (RAL). The isotopes ^'^Zn and ^^C decay by 
emitting a positron which annihilates an electron producing a pair of photons with 
an energy of 0.51 MeV. The ^^Zn and ^^C activities are determined by detecting 
the 0.51 MeV 7-rays. The ^^Zn activity shows again a quadratic power dependence 
on the laser pulse energy (Fig. 3 in^i^. 

The above mentioned experiments use different types of primary targets: solid- 
aluminum foils and deuterated clusters. The mechanisms of the ion acceleration 
from these targets are thought to be not unique, i.e., in the case of the foil targets, 
the target normal sheath acceleration (TNSA) mechanism EIIII] is responsible for 
the ion acceleration. Within this mechanism electrons on the target laser-irradiated 
surface are accelerated by the laser pulse in the direction of the laser propaga- 
tion. The electrons pass through the target, ionizing the surface contaminants, and 
form a sheath on the rear surface. The charge separation due to the motion of 
electrons causes the electrostatic field which accelerates the protons in the sur- 
face contaminants. While, in the case of deuterated cluster target irradiation, the 
deuteron-acceleration is attributed to the bound-electron expulsions from the clus- 
ters by laser-irradiations Nevertheless the obtained experimental data on the 
fusion yield show a similar quadratic dependence on the laser pulse energy. Two 
characteristic commons to those experiments are the spectra of the accelerated ions 
which have a Maxwcllian-like shape pop4 | 25 | g^j^^ ^-^^ ^^^^^ -(-j^g^-f- ^^-^c fusion reaction 
occurs in the beam-target mechanism IMIH These data have not been investigated 
comprehensively, because of the apparent difference in the ion acceleration mech- 
anisms. In this paper we analyze these data assuming hydrodynamic scaling EZIH^ 
and show that the scaling laws in fusion yields can be explained inclusively. For this 
purpose we reformulate the scaling of proton acceleration found in the maximum 
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proton energy and the conversion efficiency These two features are related 

with the temperature, i.e., the width of the Maxwellian-like spectra, of the accel- 
erated ions and the total number of the accelerated ions: both scale as functions 
of the laser power. These two parameters are, indeed, more fundamental than the 
maximum proton energy in the sense of characterizing the energy spectra of the ac- 
celerated ions and estimating the nuclear reaction rate, as is described in | Appendix] 
[B] In this context we stress the importance of determining the plasma temperature 
for the yield evaluation. 



2. Scaling of fusion yields 

We begin with the scaling features of the neutron yield from deuterium cluster 
explosions, which indeed motivated us studying the scaling of the yield from laser- 
induced nuclear-reactions. The experiments have been conducted by the group of 
the university of Texas at Austin. The neutron yield from deuterium cluster ex- 
plosions scales as a function of the laser pulse energy (E'l) the yield follows a 
quadratic power law dependence on the pulse energy. The same group has extended 
the yield measurement using different lasers with different laser parameters The 
results are summarized ^ and they are reproduced in the top panel of Fig. [l] The 
yield (iVy) shows not only the quadratic dependence but also a clear dependence 
on the pulse duration (r). By reanalyzing these data, we found that all points fall 
in a line of a quadratic power law, if {Nfr/A), where A is the laser spot size, is 
plotted as a function of the laser pulse energy. In other words, 

NfT/A^El (1) 

Dividing Eq. M by r^, the power law dependence of the yield becomes: 



Nf (El\ 



2 



This equation is suggesting that the quantity {Nj/t/A) has a power law scaling 
with an exponent of 2. The quantity {Nf/r/A) is shown in the bottom panel of 
Fig. [l] where the pulse-duration dependence seen in the top panel is disappeared. 

The scaling features found in the above experimental data of neutron yield are 
observed in the experimental data by other groups. We analyze a wider range of ex- 
perimental data by different groups worldwide together with the data by the group 
of the university of Texas at Austin. The neutron yields from three nuclear reac- 
tions D((i, n)^He, ^^B(p, n)^^C and ^^Cu(p, n)^^Zn are analyzed. Experimentally the 
data are obtained by laser-pulse irradiations on primary targets of plastic deuterium 
(CD), deuterium clusters, boron slabs and copper foils. In Tab. [l] we summarize the 
laser parameters and the observed neutron yields from the reaction D((i, n)'^He 
for various experiments. The neutron yields for the reactions ^^B(p, n)^^C and 
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Fig. 1. (Top panel) Fusion yield of the reactions D(d,n)'^He by the cluster explosions as a function 
of the laser pulse energy. The experimental data are retrieved from Ref. [19] (pulse duration 
35 fs (diamonds)) and with THOR laser at University of Texas at Austin(pulse duration 40 
fs (triangles)) [25] and with JanUSP (Calisto) laser at LLNL (pulse duration 100 fs (stars) and 1 
ps (squares)) [18]. (Bottom panel) Scaled Nf/T/A as a function of the laser power. The symbols 
for the fusion yield data are the same of the top panel. 



^^Cu(p, n)^^Zn are converted from the radioisotope activities reported EHEHl As- 
suming a beam-target fusion mechanism, the yield Nj is evaluated P^ l ^^ l as: 

Nf = {av)pTPpiAdprT (3) 

where (av) is the reaction rate per pair of particles of the reaction of interest EIIH 
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Table 1. Laser parameters, targets and the observed neutron yield for the reaction D((i, n)^He at 
various laser facilities. "The laser pulse energy is not specified in the paper and replaced by the 
ma^ximum laser pulse energy available from the laser system. ''The target thickness is not specified, 
but replaced by the absorption length in Yamanaka et al. 



Reference 




Pulse energy 


length 


Focal spot 


Target thick- 


material 


neutron yield 








(J) 


(fs) 


FWHM(/xm) 


ncss(^m) 






Pretzler et 


al 


mr 


0.2 


160. 


4.5 


200. 


C2D4 solid 


1.4x10^ 


Hilscher et 


al 


EU 


0.3 


50. 


15. 


200. 


CD2 solid 


4.xl03 


Disdier et 


al 


l34l 


30." 


300. 


5. 


400. 


CD2 solid 


10^ 


Izumi ct 


al 


EU 


50. 


500. 


5.4 


5.5 


CgDg solid 


8.8x10^ 


Belyaev et 


al 


l36l 


12. 


1500. 


15. 


200. 


CD2 solid 


10^ 


Norreys et 


al 


Ea 


20." 


1300. 


15. 


120. 


CgDg solid 


1.2x10^ 


Yamanaka et 


al 


l38l 


5. 


2 10" 


100. 


100.'' 


D2 solid 


3x10^ 








12. 


2 10" 


100. 


100.'' 


D2 solid 


5x10^ 






EH 


20. 


2 10" 


100. 


100.'' 


D2 solid 


2x10* 


Shearer et 


al 


20. 


2 10" 


80. 


100.'' 


CD2 solid 


5x10^ 








70. 


5 10" 


80. 


100.'' 


CD2 solid 


10* 


Ditmire et 


al 


l40l 


0.12 


35. 


200. 




D2 cluster 


10* 


Fritzler et 


al 


im 


62. 


1000. 


20. 




D2 cluster 


10" 


Hartke et 


al 


l42l 


0.2 


40. 


200. 




D2 cluster 


2x10^ 


Zweiback et 


al 


EU 


0.12 


35. 


200. 




D2 cluster 


5x10^ 



a and v are the reaction cross section and the velocity of the colliding nuclei, 
respectively. The term {av) takes into account that the accelerated ions have an 
energy spread which could be characterized by a temperature. We mention that {av) 
is modified from the thermonuclear reaction rate, in order to take into account that 



one of the colliding nuclei is at rest, as it is defined in Appendix B ; pT,Ppi and A 
are the number densities of the target, plasma, and the laser spot size, respectively. 
dpr is an effective target thickness. For a thin target dpr is the thickness of the 
target {d denotes the target thickness.) and for a thick target dpr is the projected 
range of the plasma ions in the target medium. In Eq. (|3]) the terms which can be 
dependent on the laser pulse energy are (erf), Ppi and dpr- If hydrodynamic flow^^H 
would be at play, we could define suitably scaled quantities which should then be 
scaling invariant depending on some dimensionless quantities such as the Reynolds 
number (Re) ESEZl Classically we can define the Reynolds number as the ratio of 
a characteristic length L and the mean free path A. We can further write L as the 
product of a typical velocity v (e.g., the thermal velocity) and a typical time r (e.g., 
the laser pulse duration). Thus 

Re = vrappi. (4) 



These are the typical quantities we have to scale, in order to obtain a universal 
behavior. We divide Eq. ([3]) by the terms which are clearly independent on the 
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laser pulse energy, i.e., (rAdpT) and get: 

Nf/r/A/d/pT = {<Jv)ppi{dpr/d) oc Re/r. (5) 

Eq. ^ has a dimension of {NumheT{N)/Time{t)/Length{L)'^/LxL^) = {N/t) and it 
is exactly the same as the scaled Eq. ^ divided by a constant (dpr)- We, therefore, 
plot n, which is defined by the l.h.s of Eq. ([5]), in Fig. [2] as a function of E^/t for 
the three nuclear reactions of interest. 

n = Nf/r/A/d/pT. (6) 

The n for all reactions shows a clear power law dependence on Ei^/t. Here we 
used d=2 mm for the cluster targets Given Eq. ^ , the power law dependence 
suggests that the reaction rate {av) for all the three reactions are similar at the 
temperature of interest. Especially for the D((i, ri)'^He reaction with the cluster 
target and with the solid target the reaction rate {av) is, in principle, identical 
for a given temperature. The power law dependence implies also that the plasma 
density {ppi) is common to all cases. In other words, the Reynolds numbers should 
be similar for all those reactions. The scattering of data observed in the plot might 
be due to some differences in the experimental set-ups not considered here, e.g., in 
the plot, we have assumed that dpr is larger than the (secondary) target thickness 
d, so that we have replaced dpr by d. Repeating the measurement under common 
experimental conditions, for instance using always very thin targets, might reduce 
somewhat the scattering of the data. In such a way one could easily scale the 
relevant features. Nevertheless the plot of Fig. [2] indicates scaling over almost as 
many as 12 orders of magnitude. Deserving special mention in Fig. [2] is the density 
of the cluster gas target which we set to 5 x 10^^ cm^'^, as given va.^^. However, 
the local density within a cluster is higher than it and it is close to the solid density 
(10^^ cm~'^)Ssl. Fig.jsjshows the H divided by the plasma critical density pcr — 10^^ 
cm^'^. In the figure the cluster target data are divided by an additional factor 8, 
to get a better coincidence between the results of the cluster target and the solid 
target. This factor could be attributed either to an ambiguity of target thickness, 
which we used 2 mm as mentioned above or to an ambiguity of the number density 
of the cluster target, or to both of them. Notice that for a real gas target the number 
of fusions would be very small since it is transparent to the laser pulses. In Fig. |3] 
the curve is the result of the best fit to all the data for the reaction D(d, 7T.)^He in 
Tab. [l] and from Ref. The curve has a power law dependence with an exponent 
2.3 and it is 

n/pcr = {o-v){ppl/pcr){dpr/d) 

- 2.1 X lO-^^EL/rf-Hcmh-^), (7) 

where E^/t is in units of 10^^ W and it is so throughout this paper unless otherwise 
specified. For the reaction ^^Cu{p,n)^^Zn the best fit is given by 

n/pcr - 7.9 X 10-i^(£;l/t)2-3(cto3s"1). (8) 
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El/t (10"W) 



63 



D(dji)"Hccl. 
D(dji)^He si. 
'Cu(pji)"zn 
B(pfl)"c 



Fig. 2. The n for the reaction D(d, n)^He with cluster targets (plusses) and with solid tar- 
gets (stars), for the reaction ^^Cu{p.n)^^Zn with copper stack (diamonds) and for the reaction 
-'^-'^B(p, 7i)^-'^C with boron slab target (squares) at laser facilities worldwide. 



with the same exponent 2.3 as the reaction D(d, n)^He. A careful observation of 
Fig. [3] reveals that the slope for the reaction ^^B(p, ri)^^C is steeper than the slope 
for the reaction 'D{d,n)^lle. The best fit for the data of the reaction ^^B(p, 7i)^^C 
gives a power law dependence with an exponent 3.5: 

U/p„ 1.3 X 10-^° {El /rf-Hcm^s-^)- (9) 

The difference in the power law dependences could be attributed to the difference 
in the reaction Q-values, as we will see in the following discussions. 



3. Origin of the scaling lavif in the laser driven fusion yield 



To investigate the origin of the slope of the scaling in the H/pcr, we reformulate 
the scaling laws for the maximum energy of laser-accelerated protons from thin 
foil targets and in the conversion efficiency of the laser pulse energy to the proton 



beam energy. Eq. (A. 4 ) derived in Appendix A gives the relation between the plasma 



temperature and the laser pulse energy 



kT ■ 



tanh 




Er 



8.7 X 10- 



0.25N 



(10) 



where me, mp,c and rj are the electron andproton masses, the speed of light and 
the conversion efficiency of the laser pulse into the hot electrons, respectively, as 
they are determined in [Appendix A| we have assumed that the plasma temperature 



July 25, 2012 0:19 WSPC/INSTRUCTION FILE hslijmpe.cite 



8 




lO" 10"^ lO"* 10-' iO"^ 10"' 10" 



El/t (10"W) 



Fig. 3. The 11 divided by the plasma critical density (pcr) for the reactions in Fig. [2] The symbols 
are the same of Fig. [2] The cluster target data are divided by a factor 8, as it is explained in the 
text. The curve is a fitting to all the data for the reaction D(Q!,n)'^He. 



is proportional to the maximum proton energy. The constant factor which should 
be multiplied to this equation to obtain the plasma temperature could be deter- 
mined by using experimental data of the energy spectra of accelerated protons. The 
spectra are measured by using either radiochromic film (RCF) stacks or a CR-39 
detector or a Faraday cup with Thomson parabola, which are placed behind the 
primary target. In Tab. [2] we summarize the laser parameters and the temperature 
of accelerated protons available from metal foil target irradiations at the VULCAN 
laser, at Calisto and Titan lasers at LLNL, at ASTRA Ti:Sapphire laser at the RAL 
and the laser facility at the Korea Atomic Energy Research Institute (KAERI). In 
Fig. |4] we plot the temperature of the plasma in the table as a function of -E^/r. 
The fitting curve shown in Fig. |4]by the solid line gives the coefficient of Eq. ( 10 ) to 



be 3.4 X 10^, i.e., the temperature of the plasma is thus related to the laser power 
(El/t) by 




kT = 3.4 X lO^^tanh^ -^J^ ^VT^^ • (H) 



Although systematical measurements are required for a more precise determination 
of this coefficient. Eq. ( 11 ) is simplified in two limits, where either the pulse duration 
(t) is much shorter/longer than the time in which a proton remains in the vicinity 
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Table 2. The number and the temperature of accelerated protons in ultrahigh-intensity laser- 
matter interactions at selected laser facilities worldwide. The number and the temperature of 
accelerated ions are from the references if those are given explicitly. Otherwise those are obtained 
by fitting the proton spectra given in the references and with an assumption that protons are 
emitted in a cone of angle 22 deg. [4]. "coated with ErHs on the back. 







r{fs) 


Target (thickness) 


N, 


kT (MeV) 


Cahsto at LLNL 




10 


100 


Au (15 ^m) 


2x10^^ 


1.4 


Titan at LLNL 


47 


150 


600 


Au (14 /im)" 


1.2x10^3 


3.3 


VULCAN at the RAL 


.20 


300 


750 


Al (10 Atm) 


IxlO^^ 


5.0 




m 


120 


1000 


Al (unknown) 


1x10^2 


2.9 


ASTRA at the RAL 


A 


0.2 


60 


Mylar (13 /xm) 


4.X10" 


0.03 






0.2 


60 


Mylar (23 ^lra) 


2.x 10"' 


0.13 


laser at the KAERI 


48 


0.3 


30 


Al (15 /^m) 


1.9x10" 


0.03 






0.3 


30 


Mylar (13 ^tm) 


1.5x10" 


0.06 




10-' 



10"^ I , , , — , , , — ' ' ' — 

0.001 0.01 0.1 1 

El/t (10" W) 



Fig. 4. Temperature of the accelerated ions from the laser-solid interaction as a function of E^/t. 
The experimental data summarized in Tab. 2 (plusses) are shown with a fit to the data with 
Eq. (A.l) (solid curve). A power law dependence of the accelerated ion temperature from the 
laser-cluster interaction [24] (dashed curve). 

of the accelerating surface charge (tq), as 

kT ^ viEL/r) [for r « 2ro) (12) 
- VvEl/t {for T » 2to). (13) 

In other words, the plasma temperature of laser-accelerated protons from thin foil 
targets depends on the laser power linearly in the limit (r ^ 2tq) or the plasma 
temperature follows square-root scaling of the laser power in the other limit (r 3> 
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2to). We note that the lasers with the pulse power higher than 100 TW correspond 
to the latter limit and that the majority of the experimental data in Fig. [3] are in 
this laser-power region. Given this condition, and that the original equation ([TT]) is 



difficult to invert, we use Eq. (13 1 to express the relation between the laser power 



and the plasma temperature. Fitting of the data gives the coefficient of Eq. (13): 

kT = T.Q^El/t. (14) 

The coefficient of this equation has a unit of (MeV/W°'^). 

In the case of laser- irradiation on deuterated cluster target, the average energies 
of emitted deuterons from D2 and CD4 cluster plasmas are known to scale as (Fig. 
4 in 1211) 

Eion - const. X E°L^'^{keV). (15) 

In the experiments the pulse length is fixed at 100 fs and the emitted ion spectra 
have been measured by using two Faraday cups, placed about 52 cm and 36 cm from 
the laser focus and at different angles. The temperature of the deuteron plasma is 
related to the ions average energy by kT — "^Eion- This leads to the power law 
dependence of the temperature of deuteron plasma on the laser power: 

kT - 1.5 X IQ-^ X (EL/rf-^^iMeV). (16) 

The curve of the power law dependence is shown by the dashed curve in Fig. |4] 
Since the plasma is most probably out of equilibrium and rapidly expanding, the 
temperature measured from the plasma ion distributions could be different from 
the temperature at the time when nuclear fusions occur, in the case of deuterated 
cluster irradiation. Nevertheless the comparison between the two curves in Fig. |4] 
shows that the power law dependence of the plasma temperature in deuterated 
cluster irradiation is clearly weaker than that of the proton plasma temperature in 
the solid target irradiation. 

Meanwhile the conversion efficiency of the laser pulse energy to the proton beam 
energy depends on the laser pulse energy linearly (Fig. 2 in^^. There the conversion 
efficiency is determined by integrating the energies of all the protons accelerated 
above 4 MeV and it is, therefore, proportional to the number of accelerated pro- 
tons (Ni): 

N, ^ (El/t). (17) 

To verify our assumption on the E^/t dependence of the number of accelerated 
ions for different pulse-length irradiations, in Fig. [s] we plot the data summarized 
m Tab. [2] as a function of E^/t. It is well known that the number of accelerated 
ions depends not only on the laser parameters but also on the density, the thickness 
and the surface treatment of the target irradiated by lasers H | 48 | 47 [ ^pj^g scattering 
of data points in the figure could be reduced by fixing the experimental conditions 
of the targets. The fitting curve to the data of the number of the laser-accelerated 
protons from thin foil targets is shown by the solid curve in Fig. [5] 

N, = 1.4 X 10^^El/t. (18) 
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Fig. 5. Number of accelerated ions from the laser-solid interaction as a function of E^/r. Exper- 
imental data at the laser facilities listed in Tab. 2 (plusses). A fitting (solid curve) to the data 
is shown and it is Ni = 1.4 X 10^^ El/t. The power law dependence on the laser power of the 
number of deuterons emitted from D2 cluster [24] (dashed curve). 



The coefficient of this equation has a unit of (W^^). 

In the case of deuterated-cluster target irradiation, the number of accelerated 
deuterons is, indeed, known to scale as a function of the laser pulse energy for 
a fixed pulse duration liSl^ i.e., Ni =4.3x lO^^J^i and Ni =3AxW^El^^, where 
the deuteron ion yield is determined using the two Faraday cups which give two 
power law dependencies. Because of the fixed pulse length (t = 100 fs) in the 
measurement, Ni is proportional to the El/t, as well. As a function of the laser 
power it is given by 

iV,; - 4. X IO^^Sl/t. (19) 

In Fig. |5] this equation is shown by the dashed curve. 

At last, the density of the plasma is, approximately, proportional to the number 
of the accelerated ions, i.e., 

Ppi ^N,^ El/t. (20) 

Now we discuss the two terms {av) and (dpr/d). An important factor that should 
be taken into account to test the scaling laws is the Q-value of the reactions. In 
the case of the three reactions of interest, the latter two reactions have negative 
Q-values, while the first one has a positive Q-value. This might change slightly 
the slopes of (crw), since the incident energy of the reaction must overcome the 
threshold energy, or alternatively the low energy ions give no contribution to the 
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yields because of the Q- value. We begin with the reaction rate at high temperatures, 
as in the case of the reaction ^^Cu(p, n)^^Zn observed at the VULCAN laser facility. 
The reaction cross section above the Coulomb barrier is approximated EIMl by 



aR = nR^^l-^, (21) 

where R, Vc and E are the sum of the nuclear radii of the colliding nuclei, the 
height of the Coulomb barrier between the colliding nuclei and the incident energy 
of the collision. Furthermore we approximate the incident energy of collision by 
the maximum proton energy Epmax, which follows the same scaling relation as the 
plasma temperature, i.e., with an intermediate value between linear and square-root 



of the laser pulse energy. (Eq.s (12) and (13)). (av) is, then, approximated by 



2E„ _ _^2.L Vc,l2x (£;^/r)0.75±0.25 



^^^"^V' E 

^ (Sl/t)"-38±0-13 (22) 

that is the leading term has a power law dependence with an exponent 0.38 ± 0.13. 
The curve for the proton projected range in a copper foil is given as a function of 
the incident energy E by the SRIM code In the energy region E > 1 MeV, it is 
approximated by 

dpr - (23) 

Here again we replace the incident energy E by the maximum proton energy. We, 
thus, obtain a relation between laser pulse energy and the projected range: 

dpr - ((i?L/r)°-^^±°-25)l-^^ ^ (^^/^)1.31±0.44^ (24) 



Finally, by using Eq.s ([20|, ([22]) and ([24|, the U/pcr for the reaction ^^Cu{p, nf^Zn 
has the iSi/r-dependence: 

n/pcr = {o-v){ppl/pcr){dpr/d) 

^ {EUrf-^'^'-'^{EL/T){EL/Ty-''^"-'' X const. 

^ (i?L/r)2-69±0-57. (25) 

The obtained (_BL/T)-dependence from this rough estimate is in accord with the 
power (El/t)'^-^, found in Eq. (|8|. Furthermore it is clearly seen that the exponent 
2.3 is closer to a limit (2.69-0.57=2.12), which corresponds to the root-square scaling 
of the plasma temperature on the laser pulse power, rather than the other limit 



(2.69+0.57=3.26). This fact supports our approximation of Eq. (11) by Eq. (13 1 



4. Formalism 

The rough estimate given at the end of the preceding section shows an excellent 
agreement with the observed power law dependence. We derive the scaling of the 
fusion yield more generally and precisely in this section. 
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4.1. Projected range and reaction rate 

In Sec. [2j the yield is expressed in terms of the projected range (dpr). That is the 
cross section is assumed to be a constant at the given incident energy, through the 
whole target. This assumption is correct, if the target is thin and the energy loss 
of the accelerated ions in the target is negligible. While for a thick target the yield 
at an incident energy (Ein) is given by integrating the beam energy down to zero 
and it is proportional to 1^33 

r° aiE') 



where Sp{E') is the stopping power of accelerated ions in the (secondary) target. 
This quantity takes into account the projectile energy loss in the target. We evaluate 
the reaction rate for the plasma ions colliding with ions at rest which is defined by 



Eq. (B.6| in Appendix B but with the substitution of (j{E) by Eq. (26): 



(avd) = I dE„, [ f P^dE'] ,/^V(i?™), (27) 



E,^ Sp{E') J \ mi 

where ijj{Ein) is the energy spectrum of the plasma ions. This gives the definition 
of (avd), which is a function of the temperature kT. In terms of (avd), Tl/pcr is 
rewritten as 

n/Pcr = {'Jv)iPpi/Pcr)idpr/d) = (avd) (ppi / Pcr)d~'^ . (28) 

The top panel in Fig. |6] shows {avd) for the three reactions as a function of the 
plasma temperature kT. As derived in the preceding section, the plasma temper- 
ature of the laser-accelerated protons from thin foil targets is related to the laser 
pulse power by a square-root dependency (Eq. ( |14P ). By using this equation we 
obtain the (i?L/T)-dependence of the (avd). In the bottom panel of Fig. [6] the 
{El I r)-dependence of the (avd) is shown for the three reactions induced by laser- 
accelerated protons from thin foil targets. Specially for the case of the D2 and 
CD4 cluster gas targets, the relation between the plasma temperature and the laser 



power is different from Eq. (14 1. For the cluster gas targets Eq. (16 1 can be used 



to obtain the (ii^L/T)-dependence of the (avd). 
4.2. Plasma density 

The density of the plasma can be dependent on the temperature of the plasma. 
To our knowledge, there is no available experimental data on the plasma tempera- 
ture (the laser pulse energy) dependence of the plasma density. Here we derive the 
laser pulse energy dependence of the plasma density in two limits, using the power 
law dependence of the number of accelerated ions (Eq. ([T7|). The first limit is the 
case where the volume of the plasma (Vpi) has no £^i/r-dependence. The plasma 
density is related to the laser pulse power linearly: 

Pp/ = ^ ^ El/t. (29) 
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D(d^) He X 
"b(pji)"c □ 
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D(d,n} He x 
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"Cu(p.n)"Zn 
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10" 10' 



Fig. 6. (Top panel) The (avd) defined by Eq. (27) as a function of the plasma temperature. For 
the reactions D(d, n)^He (crosses), ^^3{p,n)^^C (squares) and ^•^Cu(p, n)^'^Zn (triangles). The 
cross sections are retrieved from the NACRE [44], EXFOR [55] databases. The stopping power 
data are obtained using the SRIM code [51]. (Bottom panel) The same with the top panel but as 
a function of the laser power, derived with Eq. (14). 



The other limit is a case where the volume of the generated plasma has the 
laser-pulse power dependence, e.g., Vpi is expressed as a product of the laser 
focal spot size and the range covered by the accelerated ions in a period of 
the laser pulse duration, i.e., Vpi = Ad. This range can be approximated by 
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cT = 2Epmax /iTipT. Given the laser-pulse power dependence of the max- 



imum proton energy Epmax [El/t) (Eq.(22)), the volume of the plasma is 
approximated by Vpi = AT{EL/TY'^'^^ . The approximation is justified because the 
energy loss of the energetic ions in the plasma is known to be negligible Substi- 
tuting this, the plasma density is expected to depend on the laser pulse power: 

_ {El It) 



Ppi 



(30) 



At{El/t)0-2^' 

Which one of the two limits is closer to the real case might be inferred from the 
scaled Il/pcr- The scaled H/pcr in Eq.s Q and ^ is rewritten in terms of 

{(Tvd) as 

^/Pcr ^ {tyvd){ppi/pcr)d^^- (31) 



In the first limit (Eq. (29l), this equation becomes 

n/p,r = {(Jvd)iEL/T){pcrd)-\ (32) 

where the last term of the right-hand-side is a constant. The iJ^/T-dependence 



comes from the first two terms. In the second limit, (Eq. (30l), Eq. (31 1 gives: 

{avd){EL/Tf-'^p,rAdTy\ 



(33) 



We evaluate Eq. (32) and Eq. (33 



laser power dependence of Eq. 



numerically. It should be emphasized that the 



32 1 and Eq. (33 1 are derived using Eq. (13) 



4.3. Discussions 

In the top panel of Fig. [t] the /r-dependent term of Eq. ( p2| ) is shown by 
bigger symbols connected with thick lines for the reactions D{d,n)^lle (crosses), 
^^B(p, n)^^C (squares) and ^■^Cu(p, n)^'^Zn (triangles). The /'''-dependent term 
of Eq. ( [SS] ) is shown by smaller symbols connected with thin lines for the same 
reactions. The thin straight dashed line represents the curve of (El/t)^-^, i.e., the 



observed scaling relation of experimental data for the reaction D{d, n)'^He (Eq. ([7|) 
For the reaction D(d, n)'^He, it is clearly seen that the curve with the exponent 2.3 



is close to the estimate by Eq. (32), while the estimate by Eq. (33) gives an expo- 



nent less than 2. The comparison with the curve with the exponent 2.3 suggests 
that the plasma density depends linearly on the laser power, i.e., the volume of 
the generated plasma has no i?L/T-dependence. The same apphes for the reaction 
^^Cu{p,nf^Zn in the laser power region 3 x lO^^ <^ El/t < 3 x lO" W, where 
one finds the experimental data. While for the reaction ^^B{p,n)^^C, in the laser 
power region 3 x 10^^ W < El/t < 5x 10^"^ W, the curve given by Eq. (32 ) has the 
(i?i/r)-dependence of {El/t)'^-^. This is in good agreement with the laser power 
dependence of the exponent 3.5 found in Eq. ([£]). In the bottom panel of Fig. [7| 
the thick curves given by Eq. ( |32[ ) for the three reactions are compared with the 
experimental data of 11 in Fig. l2P[n the figure Eq. ( 32 ) is multiplied by a constant 



to be compared with each of their experimental data. The laser power dependences 
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Fig. 7. (Top panel) The laser power dependence of H/pcr evaluated in two limits for the reactions 
D((i,n)'^He (crosses), ^^B(p, n)^^C (squares) and ^'^Cu(p, n)^'^Zn (triangles). One limit is calcu- 
lated by Eq. (32) (bigger symbols) and the other limit is given by Eq. (34) (smaller symbols). The 
curve of (E^/r)^'^ is given by the straight dashed line. The abscissa is limited in the same region 
of Fig. 3. (Bottom panel) The thick curves in the top panel are compared with the experimental 
data of the three reactions induced by laser-solid target interaction. 11 of the experimental data 
is shown by points as they are defined in Fig. 2. 



of {avd)Ni derived under an assumption of the linear dependence of the plasma 
density on the laser power (Eq. (32|) for the three reactions are in accord with 11 
of the experimental data. The steep rises of the curves given by Eq. (32) for the 
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10^ 




10"'' 10"-' 10^ 10"-' I0-' 10-' 10^ 

El/t (10"W) 



Fig. 8. Same with the bottom panel of Fig. [t] but for the reaction D(ii,n)''He induced by the 
D2 cluster gas target irradiatio n. T he laser power dependence of {avd)Ni derived with kT ~ 
0.015 X (El/t)°'^^ (MeV) (Eq. |l6|) for the reaction B{d,n)'^Ke (solid line) with the D2 cluster 
gas target. The experimental data of 11 are shown by pluses. 



reactions ^^B{p,n)^^C and ^'^Cu(p, n)^'^Zn in the laser power region less than 5 x 
IQi^ \Y are due to the negative Q- values of the reactions. 

In the case of laser-irradiation on the deuterated cluster target, the power law 
dependence of the deuteron plasma temperature on the laser power is expressed 
by Eq. 
in Fig. 



We evaluate Eq. (32) numerically by using Eq. (16) and it is shown 



by the solid curve together with 11 of the experimental data (shown by 
pluses). The laser power dependences of {avd)Ni derived by assuming the linear 
dependence of the plasma density on the laser power (Eq. (32 )) is again in excellent 
agreement with 11 of the experimental data using deuterated cluster targets. 

It deserves special mention that the derivation of the scaling relation relies on 
the fact that energy spectra of accelerated ions are characterized by a Maxwellian- 
like shape with a temperature. One should verify the shape of the spectra of the 
accelerated ions, before extending the result of the scaling law to the laser power 
beyond PW (El/t > IQi^ W). 



5. Summary 

We have derived scaling laws in the fusion yields resulting from intense laser irra- 
diations on either solid or cluster gas target, as a function of laser parameters. The 
origin of the power scaling has been studied and it is attributed to the laser power 
dependence of three terms: the reaction rate, the density of the plasma and the 
projected range of the plasma particle in the medium. The density of the plasma 
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is derived to scale linearly as a function of the laser power. The resulting scaling 
relations can be used, e.g., to estimate the yield of positron emitters produced by 
deuteron-induced reactions by means of laser-accelerated deuterons. The obtained 
scaling law enables estimating the fusion yield for a nuclear reaction which has not 
been investigated by means of the laser accelerated ion beams. 



Appendix A. Maximum proton energy and plasma temperature 

The maximum proton energy of laser-accelerated protons from thin foil targets is 
known to scale as a function of the laser power 

Epmax = Eoo tanh^ (2^) ' ^^''^'^ 

where £'00 and tq are the potential barrier in which the hot electrons are confined 
and the time in which a proton remains in the vicinity of the accelerating surface 
charge; t is the pulse duration as before. E^o and tq can be written as a function 
of the laser power: 



where me, c and rj are the electron mass, the speed of light and the conversion 
efficiency of the laser pulse into the hot electrons, respectively, and 



where rup is the proton mass. Substituting the equations (A.2) and (A. 3 1 into 



Eq. (A.l), we obtain the laser pulse energy dependence of the maximum proton 
energy. The isothermal fiuid model gives the relation between the maximum proton 
energy and the hot electron temperature EMU Taking into account this relation, 
we assume that the plasma temperature is proportional to the maximum proton 
energy. 




Appendix B. Reaction rate per pair of particles 

The thermonuclear reaction rate is evaluated byl^ 

(av) = / a{v)v(f){v)dv^ , (B.l) 



where (piv) is the relative velocity spectrum of a pair of ions and it is given by a 
Maxwellian-distribution at the temperature kT: 
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denoting the reduced mass of the colliding ions in MeV as fi. Here we write the 
reaction cross section as a function of the velocity to keep the consistency in the 
velocity integral, instead of as a function of the incident energy. Since we consider 
the collisions between the ions in the plasma and the nuclei in the target (HT) , one 
should take into account that one of the colliding ions is at rest in the laboratory 



frame. Therefore the velocity spectrum Eq.(B.2) is modified as: 

,2 " 



where mi is the mass of the ions in plasma, instead of the reduced mass. One can 
define an effective temperature as, 

k^eff ^ JLkT. (B.4) 

mi 



Using this, the velocity spectrum Eq. (B.3 1 is converted to the energy spectrum. 
The energy distribution is therefore: 

The reaction rate for the plasma ions colliding with ions at rest is calculated by 

lav) ^ ! cT{E)J—ME)dE. (B.6) 
J \ mi 

In 

Fig. [ojthe reaction rates are shown for D((i, ?i)'^He (crosses), ^^B(p, rt)^^C (squares) 
and ^■^Cu(p,n)^'^Zn (triangles), respectively. The thermonuclear reaction rates for 
the three reactions are given by the small symbols. The data are retrieved from 
the NACRE-compilation ^ for the first two reactions and from a database of the 
Hauser-Feshbach statistical model calculations 1^ for the reaction ^■^Cu(p,n)^'^Zn. 
While the reaction rates between an ion in the plasma and an ion at rest are shown 
by bigger symbols. The quantity of interest is the latter and the reaction rates in 
plasmas (smaller symbols) are shown just for comparison. The reaction rates for 
the reactions ^^B(p, n)^^C and ^^Cu(p, 7i)^'^Zn show steep rises at the temperature 
around 100 kcV. This is due to the negative Q-values for the two reactions. 

Appendix C. Effective Temperature 



In Eq. (B.4 1, the factor in front of the temperature is the same factor which connects 
the relative energy of two colliding particles (E) in the center-of-mass system and 
the relative energy in the laboratory system (Eiab)- 

E = —Eiab, (C.l) 
mi 

where Eiab — mif^/2. It is the same either we work out in the laboratory system or 
in the center-of-mass system. We take the center-of-mass system and multiply the 
temperature by the conversion factor of the energies, for convenience, so as that the 
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Fig. 9. Reaction rate per pair of particles for D(d, n)^He (crosses), ^^B(p,n)^^C (squares) and 
®^Cu(p, n)®^Zn (triangles). The smaller symbols are for the thermonuclear reaction rates [44,52]. 
The bigger symbols are for the reaction rate between an ion in the plasma at the specified tem- 
perature and an ion at rest. 



energy E in Eq. (B.5) is still the energy in the center-of-mass system. By using the 
effective temperature determined by Eq. (B.4|, the velocity distribution (Eq. (B.3)) 
is 



(j)HT{v)dv'^ 



( A* 



V 27rfcT, 



exp 



2kT, 



Anv'^dv. 



(C.2) 



This is exactly the same with the velocity distribution of Eq. (B.2 1 but the temper- 
ature is replaced by the effective temperature. The energy distribution Eq. (B.5) is 
obtained by substitution of 

M 2 



E ■ 



(C.3) 



in Eq. (C.2 1. 
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